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Abstract 



In this paper, we present results of numerical lattice simulations of two- 
flavor QED in three space-time dimensions. First, we provide evidence that 
chiral symmetry is spontaneously broken in the chiral and continuum limit. 
Next we discuss the role of an external magnetic field B on the dynamically 
generated fermion mass. We investigate the 5-dependence of the condensate 
through calculations with dynamical fermions using the non-compact formu- 
lation of the gauge field, and compare the results with those of a comparable 
study using the quenched approximation. 
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j£ ; 1 Introduction 

The intriguing phenomenon of high T c superconductivity has inspired several 
models, consisting mostly of elaborations of a doped version of the Hubbard 
model. The model presented in |] claims the generation of superconduc- 
tivity through a "mixed Chern-Simons term" without the need to introduce 
a local gauge symmetry breaking condensate. The model invokes the idea 
of electron substructure (spin-charge separation) considering the low-energy 
excitations of the system to be bosonic degrees of freedom called spinons and 
fermionic ones called holons. The role of the phonon interactions crucial to 
the formation of Cooper Pairs in BCS superconductors is now played by a 
"statistical" gauge field, which arises from the underlying dynamics of the 
holons hopping between lattice sites. More specifically, the model is based 
on a so-called T3-U(l) interaction of the fermions with the gauge field. 

An alternative form of the model has been proposed in ||. The model 
is meant to represent the naive continuum limit of the effective low-energy 
Hamiltonian of the \axge-U, doped Hubbard model. The action of the model 
reads: 

C = -\{F, v f - ^tr(^) 2 + VDtffi (1) 



where _D M = <9 M - ig\A s - ig2"T a W a ^ - i- c A e * 1 , with A$ the statistical U(l) 5 
gauge field, is the gauge potential of a local SU(2) group, and the 
corresponding field strength. The \I/ fields are two-component spinors. The 
difference from the previous model lies in the fact that now the full gauge 
symmetry is U(l)s ® SU(2), where we have denoted by U(l)s the abelian 
part of the symmetry. We call it a "statistical" gauge field since from the 
physical point of view it describes statistical properties of the model, such as 
the doping, and not real electromagnetism. 

The strongly-coupled U(l)s interaction promotes the spontaneous gener- 
ation of a fermion - anti-fermion condensate resulting in a mass gap. Argu- 
ments based on energetics prohibit the generation of a parity-violating but 
SU(2) gauge invariant term ||, in favour of a parity-conserving condensate, 
which necessarily breaks the SU(2) group down to a r 3 -U(l) sector Q. The 
mass gap thus breaks the SU(2) part down to an abelian sub-group which is 
the analog of the r 3 -U(l) in the model of The symmetry breaking pattern 
here follows the scheme: U(l) 5 <g> SU(2) -> U(l) s <g> U(1) T3 . 

In equation ([!]) we have also allowed for an additional U(l) em coupling, 
which describes the coupling to an external electromagnetic field A e ^ 1 . Actual 
superconductivity is obtained upon coupling the system to such external elec- 
tromagnetic potentials. We now discuss the superconducting consequences of 
the above dynamical breaking patterns of the SU(2) group. Upon the opening 
of a mass gap in the fermion (hole) spectrum, one obtains a non-trivial result 
for the following Feynman matrix element: S a = (W^\J u \0) , a = 1, 2, 3, with 
J M = \j/7 M \I/ the fermion-number current. Due to the colour-group structure, 
only the massless W 3 gauge boson of the SU(2) group, corresponding to the 
T3 generator in two-component notation, contributes to the matrix element. 
The non-trivial result for the matrix element S 3 arises from an anomalous 
one-loop graph, and is given by 0. [lj : 

S 3 = (W 3 \J„\0) = ( S gnM) W -^ (2) 

where M is the parity-conserving fermion mass generated dynamically by the 
U(l)5 group. As with other Adler-Bell-Jackiw anomalous graphs in gauge 
theories, the one-loop result (fj) is exact and receives no contributions from 
higher loops [|]]. This unconventional symmetry breaking (Q), does not have 
a local order parameter j|, [I]] since its formation is inhibited by strong phase 
fluctuations, thereby resembling the Kosterlitz-Thouless mode of symmetry 
breaking. The massless gauge boson W 3 of the unbroken r 3 -U(l) subgroup 
of SU(2) is responsible for the appearance of a massless pole in the electric 
current- current correlator, which is the characteristic feature of any super- 
conducting theory. As discussed in ref. |[], all the standard properties of a 
superconductor, such as the Meissner effect, infinite conductivity, flux quan- 
tization, the London limit of the action etc. are recovered in such The 
field W 3 , or rather its dual (ft defined by <9 M = e^ p d u W 3 , can be identified 
with the Goldstone boson of the broken U(l) em electromagnetic symmetry 
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If an external magnetic field is coupled to the charged excitations about 
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the superconducting state then new phenomena appear. In refs. [0, Q it 
is argued that an external magnetic field in this (2+l)-dimensional system 
induces the opening of a second superconducting gap at the nodes of the 
ci-wave gap, in agreement with recent experimental results on the behaviour 
of the thermal conductivity of high-temperature cuprates under the influence 
of strong external magnetic fields ||. 

In this work we will study the phenomenon of the enhancement of the 
fermion condensate in the presence of an external magnetic field. In this 
respect it is not important to gauge the SU(2) group. We study QED 3 , a 
model which in the continuum has a local U(l) symmetry and a global SU(2) 
symmetry. However, the lattice discretization we employ in order to perform 
non-perurbative calculations restricts the global symmetry to a U(l) "chiral" 
symmetry, which in turn is totally broken by the fermion - anti-fermion 
condensate The question of whether this symmetry is spontaneously 

broken in QED 3 , ie. whether (^vp) 7^ in the chiral limit bare fermion 
mass m — > and the continuum (weak coupling) limit g — > 0, has been 
an important issue in non-perturbative field theory for many years. Initial 
studies based on Schwinger-Dyson (SD) studies using the photon propagator 
derived from the leading order 1/Nf expansion, where Nf is the number of 
fermion flavors, suggested that for Nf less than some critical value Nf c the 
answer is positive 0, with Nf c ~ 3.2. The model in the limit Nf — > Nf c 
is supposed to undergo an infinite-order phase transition [||. Other studies 
taking non-trivial vertex corrections into account predicted chiral symmetry 
breaking for arbitrary Nf More recent studies which treat the vertex 
consistently in both numerator and denominator of the SD equations have 
found Nf c < 00, with a value either in agreement with the original study 
[|IIJ, or slightly higher Nf c ~ 4.3 [JL1]J . Finally a recent argument based on a 
thermodynamic inequality has predicted Nf c < | [|i~2 |. 

There have also been numerical attempts to resolve the issue via simu- 
lations of non-compact lattice QED3. Once again, opinions have divided on 
whether Nf c is finite and ~ 3 |13|| , or whether chiral symmetry is broken for all 
Nf U|] . The principal obstruction to a definitive answer has been large finite 
volume effects resulting from the presence of a massless photon in the spec- 
trum, which prevent a reliable extrapolation to the thermodynamic limit. 
In this context it is also worth mentioning the three-dimensional Thirring 
model, which has a current-current interaction mediated in the UV limit by 
a propagator of identical form in the large- iVj limit to that of QED 3 ; more- 
over chiral symmetry breaking for small Nf is predicted by the same set of SD 
equations as those of QED3 ||15| . This time, however, the critical behaviour 
is governed by a UV-stable fixed point rather than an IR one, and numerical 
studies are more controlable, yielding values of Nf c ~ 5 [16]. 



The rest of the paper is organised as follows. In section || we outline the 
formulation of continuum QED 3 in Euclidean metric, and identify its global 
symmetries. In section |3| we do the same for the lattice version of the model. 
We also review the formulation of a homogeneous background magnetic field 
B in abelian lattice gauge theory. Our numerical results are presented in 
section [|. First for B = 0, we attempt an extrapolation to both continuum 
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and chiral limits, and observe evidence for a non-vanishing (though unex- 
pectedly small) dimensionless chiral condensate, suggesting Nf c > 2. We 
are also able to estimate the combination of lattice sizes and bare masses 
required for quantitative accuracy in further work. This work is described 
in section |4j]. In section [4.2| we turn to non-zero magnetic field B, and 
study the resulting enhancement of the condensate. This is the first lattice 
simulation of QED3 with B 7^ using dynamical fermions; the problem has 
previously been studied in the quenched approximation JF/J. In particular 
we compare the behaviours of quenched versus dynamical fermions and cal- 
culate the condensate as a function of the magnetic field, and find different 
behaviours both in the strong and the weak field regimes, and in strong and 
weak U(l)s coupling regimes. We conclude with a study of the dependence 
of the condensate in the strength of the statistical gauge interaction in the 
chiral limit, and uncover an important finite volume effect. 



2 The model in the continuum and its sym- 
metries 

The three-dimensional continuum Lagrangian describing QED 3 is given (in 
Euclidean metric, which we use hereafter) by 



£ = -F^ + ViD^i + m%% (3) 

where = — igA^ — ieA^, and F^, are the corresponding field 
strengths for an abelian U(l)s statistical gauge field A^ and a background 
U(l) em gauge field A^, respectively. The fermions % = l,...,Nf, are 
now four-component spinors. We note that may be written in terms of 
the two-component spinors used in section |1] as 



»« - I' 1 ■ w 




A convenient representation for the 7^, fi = 0, 1, 2, matrices is the reducible 
4x4 representation of the Dirac algebra in three dimensions (in Euclidean 
metric we choose hermitian Dirac matrices): 

where <Ji are 2x2 Pauli matrices. Then the Lagrangian (|3]) decomposes into 
two parts, one for and one for \l/j2, which will be called "fermion species" 
in the sequel. The presence of an even number of fermion species allows us to 
define chiral symmetry and parity in three dimensions 0, which we discuss 
below. The bare mass m term is parity conserving and has been added by 
hand in the Lagrangian (^). This term is generated dynamically anyway via 
the formation of the fermion condensate (\E , \I>) by the strong U(l)s coupling. 
Here we include such a term explicitly, since it is necessary from the technical 
point of view. 
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As is well known ]7| there exist two 4x4 matrices which anticommute 
with fi = 0, 1, 2: 

73 = (" n ) > 75 = i ( n ) , (6) 







■( 
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J), 75 = 


7-i 


J) 


are 


2x2 matrices. 


These 


are 



'chiral' symmetry for the massless fermion theory: 

* -> exp(z6 l 7 3 )\l/, * — > ^exp(^7 3 ); 

* -> exp(zcJ7 5 )^, * -> ^exp(iu;75). (7) 

Note that these transformations do not exist in the two-component repre- 
sentation of the three-dimensional Dirac algebra, and therefore the above 
symmetry is valid for theories with even numbers of fermion species only. 

For later use we note that the Dirac algebra in (2 + 1) dimensions satisfies 
the identities: 



l o 

-1 



rf^'f'f = -5 pfT 7° - 3r 3 e pCT 
YYl a Y = -35 pa - r 3 7°e p<T 

^YYYi 11 = -5 P<7 Y - 5 pT Y + 5 aT Y (8) 

which are specific to three dimensions. Here the Greek letters are spacetime 
indices, and repetition implies summation. 

Parity in this formalism is defined as the transformation: 

P : ^(x ,^ 1 ,^ 2 ) -> -ry 3 7i^(x°, -x 1 ,^ 2 ); 

^(xV 1 ,^) -> ^(x°,-a; 1 ,a; 2 )i7i7 3 , (9) 

and it is easy to see that the parity-invariant mass term corresponds to 
the species having masses with opposite signs [J/J, whilst a parity- violating 
one corresponds to masses of the same sign. 

The set of generators 

£ = {1, 73 ,75,t 3 } (10) 

form [18, U a global U(2) ~ SU(2)(g)U(l) symmetry. To see this define 
\1/ = \I/t 3 . The U(2) symmetry is then 

V->UV -> W f C/ G U(2). (11) 

The identity matrix 1 generates the U(l) subgroup, while the other three 
form the SU(2) part of the group. The currents corresponding to the SU(2) 
transformations are: 

jj = %r* r = 73 , 75 ,r 3 (12) 
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and are conserved in the absence of a fermionic mass term. It can be readily 
verified that the corresponding charges Qr = / <i 2 x\I/7or\l/ lead to an SU(2) 



algebra [18 



[QsiQs] = 2iQ T3 , [Q 5 ,Qr 3 ] = 2lQ 3 , [Q 7 *,Q 3 ]=2iQ S . (13) 

In the presence of a mass term, these currents are not all conserved: 

d^Jl = 2m^r^, (14) 

for T = {73, 75}, while the currents corresponding to {1, r 3 } are always con- 
served, even in the presence of a fermion mass. The situation parallels the 
treatment of the SU(2)®SU(2) chiral symmetry breaking in low-energy QCD 
and the partial conservation of axial current (PCAC). The bilinears 

At = ^73^, A 2 = #75^, A3 = 

Bin = *7„73*, B 2li = ^7/,75^, B 3fM = ¥7^3*, fi = 0, 1, 2 (15) 

transform as triplets under SU(2). The SU(2) singlets are 

A4 = ^r 3 ^, = V-yfl (16) 

i.e. the singlets are the parity-violating mass term, and the four-component 
fermion current. 

We now notice that in the case where the fermion condensate A3 is gen- 
erated dynamically, energetics prohibit the generation of a parity-violating 
SU (2) -invariant term [[|, and so a parity-conserving mass term necessarily 
breaks § the C/ S (1)®SU(2) group down to a f/ 5 (l)®U(l) T3 group 0. 

Finally, we note that were the global symmetries described here to be 
promoted to local ones, as they are in both the r 3 -U(l) model of and 
the C/s(l)®SU(2) model of 0, then the measure of the Euclidean functional 
integral (i.e. the fermion determinant) would not be positive definite. Pos- 
itivity of the theory defined by @ follows because there exist matrices 73, 
7 5 which anticommute with D^^; hence the imaginary eigenvalues of this 
antihermitian operator necessarily occur in conjugate pairs. This property 
no longer holds once, e.g., the U(1) T3 symmetry is gauged. Thus this class 
of models suffer from a 'sign problem' in Euclidean metric, a feature shared 
with other potentially superconducting systems such as QCD at high baryon 
number density [|19l . 



3 Lattice formulation 

We now proceed with a description of the lattice formulation of the problem. 
The lattice action using staggered lattice fermion fields \i X is given by the 
formulae below: 

5 = ^ E F^(x)F^(x) + £ X (x)Q(x, x') X (x') (17) 

x,n<u x ,x' 
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F^(x) = a$(x) + a s v {x + fi)- a£(x + v) - a s v {x) 
Q(x,x') = m5 X)X i + -^■n ll {x)[8 x >p+ tL U xv y xil - (W-a^-a^K-aJ- 

/ A 4 

The indices x, x' consist of three integers (xi, X2, £3) labelling the lattice 
sites, where the timelike direction is now considered as the third. Since the 
gauge action F 2 is unbounded from above, fllTD defines the non-compact for- 
mulation of lattice QED. The 77^ (x) are Kawamoto-Smit phases (— l) :ClH ^m-i 
designed to ensure relativistic covariance in the continuum limit. For the 
fermion fields antiperiodic boundary conditions are used in the timelike di- 
rection and periodic boundary conditions in the spatial directions. The phase 
factors in the fermion bilinear are defined by U xli = exp(iaf ), V xli = 
exp(ia^), where represents the statistical gauge potential and 
the external electromagnetic potential. In terms of continuum quantities, 
a xfl = agA^(x), the coupling (3g = -5^, and a^ft = aeA e ^ ct (x), where a is the 
physical lattice spacing. The coupling e to the external electromagnetic field 
is dimensionless. 

The numerical results in this paper were obtained by simulating the action 
fllTD using a standard Hybrid Monte Carlo (HMC) algorithm. The form of 
( |T7D permits an even-odd partitioning, which means that a single flavor of 
staggered fermion can be simulated. In the long wavelength limit, this can 
be shown to correspond to Nf = 2 flavors of the four- component fermions 
\]/ considered in section [| [p0|| . Chiral condensation (tyty) ^ would then 
result in a pattern U(2iVy) — > XJ(Nf) ®U(iVy) T3 . For non-zero lattice spacing 



a, however, the global symmetries for m = are only partially realised flic 
In this case the symmetry is 

Xo y Xo& Xe y 6 Xe 

Xe^Xee 113 Xo-^e^Xo, (18) 

where Xo/e denotes the field on odd (i.e. e(x) = (— l^i+^+^a — _1) anc [ 
even sublattices respectively. Chiral symmetry breaking characterised by 
(xx) ^ now has the pattern U(iV) <g> U(iV) -> U(iV) for N flavors of 
staggered fermion. 

We should now construct a lattice version of the homogeneous magnetic 



field. This has already been done before in |2T| in connection with the abelian 
Higgs model. We follow a slightly different prescription, which we describe 
below H, and which is closely related to a prescription for introducing a 



uniform topological charge density in two dimensions p2 . 

Since we would like to impose an external homogeneous magnetic field 
B = V x a ext in the (missing) z direction, we choose the external gauge 
potential in such a way that the plaquettes in the X1X2 plane equal B, while 
all other plaquettes equal zero. One way in which this can be achieved is 
through the choice: a^ xt (xi,x 2 ,x 3 ) = 0,for all xi, x 2 , x 3 , and 

a, [x 1 ,x 2 ,x 3 )-^_s {N + 1){x2 _ 1 ^ Xi = N: 

n ext (r r \-\ +1(^-1), x 2 ^N; 

a 2 (^2,23)- j + | (jV + 1)(xi _ 1); X2 = N _ (19) 
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The lattice extent is assumed to be N x iV in the 21X2 plane. It is trivial 
to check that all plaquettes starting at (xi, £2, 23), with the exception of the 
one starting at (N, N,Xs), equal B. The latter plaquette equals (1 — N 2 )B = 
B — (N 2 B). One may say that the flux is homogeneous over the entire £1X2 
cross section of the lattice and equals B. The additional flux of — (N 2 B) can 
be understood by the fact that the lattice is a torus, that is a closed surface, 
and the Maxwell equation V • B = implies that the total magnetic flux 
through the X1X2 plane should vanish. This means that, if periodic boundary 
conditions are used for the gauge field, the total flux of any configuration 
should be zero, so the (positive, say) flux B, penetrating the majority of the 
plaquettes will be accompanied by a compensating negative flux — (N 2 B) 
located on a single plaquette. This compensating flux should be "invisible" , 
that is it should have no observable physical effects. This is the case if the 
flux is an integer multiple of 2n : N 2 B = n2n =^> B = rij^, where n is an 
integer. Thus we may say (disregarding the invisible flux) that the magnetic 
field is homogeneous over the entire cross section of the lattice.^] The integer 
n may be chosen to lie in the interval [0, ^-], since the model with integers 
between and N 2 is equivalent to that with integers taking on the values 
N 2 — n. It follows that the magnetic field strength B in lattice units lies 
between and ir. 

The question of how to identify the continuum magnetic field strength 
Bphys requires careful discussion. In the physical problem we wish to address 
the magnetic field is actually the z component of a four-dimensional vector 
field, implying a relation B = ea|V x A ext \ = ea 2 B p h ys , where e is the 
electronic charge and is dimensionless. Therefore the combination (3qB is 
dimensionless. The physical field may go to infinity letting the lattice spacing 
a go to zero, while B is kept constant. 

On the lattice the magnetic field is determined by the number n, but the 
allowed values of this number depend on the lattice size N. On the contrary, 
the parameter B does not depend much on the details of the discretization, 
so it is natural to choose it as the most natural lattice unit for the magnetic 
field. In figure [I] we show the fermion condensate versus n. The solid squares 
represent the condensate for a 12 3 lattice for several values of n. We observe 
that the condensate increases substantially with n, since different n yield 
different values for B. The triangles represent data produced by taking each 
n with an appropriate volume, such that B remains almost constant; the 
resulting (n, N 3 ) pairs turn out to be: (3,8 3 ), (7, 12 3 ), (12, 16 3 ), (19, 20 3 ). 
It is obvious that in this way we get essentially the same condensate for all 
n. 

An important remark is that the magnetic field should not be allowed to 
grow too big in lattice units, since then the perturbative expansion of the 
expansions exp(za^f) would yield significant B 2 , B 3 , . . . contributions with 
the accompanying vertices, in addition to the desired terms which are linear 
in B. A trivial estimate of the critical field strength is obtained from the 

lr To check this translational invariance we measured the fermion condensate at every 
point in the x±X2 plane. The results were the same at all points within the error bars, 
confirming homogeneity. 
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Figure 1: Fermion condensate versus n for m — 0.1 and (3g = 0.30. The solid 
squares are the results for a 12 3 lattice. The triangles are data for which the 
volumes are different for each n. 

demand that the cyclotron radius corresponding to a given magnetic field 
should not be less than (say) two lattice spacings. This calculation yields 
B < ~. Of course the above limitations apply strictly only to the case where 
the statistical gauge field has been turned off; in the interacting case, one 
does not really know whether there exists a critical magnetic field after which 
discretization effects are important. With this remark in mind, we depict in 
the figures of the following sections the results for the whole range of the 
magnetic field, from to n. 



4 Simulation results 

4.1 B = 

It is instructive to begin by studying the case B = 0. We have performed 
simulations on 12 3 and 16 3 lattices for a range of gauge couplings (5q G 
[0.5, 0.8] and bare fermion masses m G [0.01, 0.05]. For each parameter set we 
accumulated 0(2000) HMC trajectories of mean length 0.9, with timesteps 
ranging from dr = 0.05 on 12 3 with m = 0.05 to dr = 0.025 on 16 3 with m = 
0.01. The chiral condensate (xx) was estimated by stochastic means using 10 
noise vectors every two trajectories. Our data is shown in Fig. It is clear 
that finite volume effects are significant; this can be attributed to the presence 
of a massless particle, the photon, in the spectrum of the theory. Large finite 
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volume effects are the main source of difficulty in numerical studies of QED 3 
[ |T3| , (HTj . We will not attempt to extrapolate to the thermodynamic limit, 



but instead focus on the 16 3 results. 
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Figure 2: 0q{xx) vs. Pg on 12 3 (filled) and 16 3 (open) lattices. 
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Figure 3: /3 G (xx) vs. Pc m from a 16 3 lattice. 
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Note that in Fig. |2| we have chosen to plot the dimensionless combination 
Pg(xx)- Since QED 3 is an asymptotically- free theory (recall g has positive 
mass dimension), its continuum limit lies in the limit (3g — > oo. To see 
whether lattice data taken with a ^ are characteristic of the continuum 
limit, it is helpful to plot dimensionless quantities against each other; results 
from different (3q should fall on a universal curve if continuum physics is being 
probed. We show a plot of (3q(xx) vs - Pcm m Fig- 0- The data do indeed fall 
approximately on a straight line, the main departures being at the strongest 
coupling (3q = 0.5, and the smallest mass m = 0.01, which might plausibly 
be attributed respectively to lattice artifacts and finite volume effects. If we 
exclude these points, a straight line fit of the form 0q{xx) = 0.00327(15) + 
1.454(6)/3(3m is achieved (if the equivalent points from 12 3 are used, the fitted 
line is 0.00140(16) + lA3Q(3 G m). 

The main result of these fits is that the value of /3 G (xx) extrapolated to 
the chiral limit is unusually small, but significantly different from zero. If 
the linear form really does characterise the continuum limit, then this result 
is robust even on our admittedly small volumes. The implication is that 
chiral symmetry is broken in the the IR limit for Nf = 2 flavors of fermion 
in QED3, and hence Nf c > 2. This is consistent with the findings of most 
studies using the gap equation, but contradicts a recent bound Nf c < | 
derived by counting massless degrees of freedom via the absolute value of the 



thermodynamic free energy / in both IR and UV limits [12|, and demanding 



fiR < fuv- (20) 



A possible resolution follows from the discussion below eqn. (|THp ; the global 
symmetry of the lattice model is smaller than that of continuum QED 3 , 
resulting in a smaller number of Goldstone bosons in the chirally broken 
phase. For N flavors of staggered fermion, there are 8N = 4Nf fermionic 
degrees of freedom in the UV limit, and iV 2 Goldstone boson degrees of 
freedom in the IR. Following the argument of |I2[ , we arrive at the bound 
N c < 6 for the lattice model, which clearly is not threatened by our results. 
To substantiate this it will be necessary to show that the strongly-coupled 
dynamics of the IR limit is governed by the symmetries of the lattice model 
rather than the continuum one |IT6f . 



Next we discuss longitudinal and transverse susceptibilities, defined re- 
spectively as the integrated propagators in scalar and pseudoscalar meson 
channels: 

xi = J2(xx( Q )xx(%)}; Xt = J2(xzx(fyx£x(x)} ■ (21) 

X X 

The longitudinal susceptibility xi has contributions from diagrams with both 



connected and disconnected fermion lines |Tj|; we find that the connected 
contribution is overwhelmingly dominant throughout our parameter space. 
The transverse susceptibility is most conveniently calculated via the axial 
Ward identity Xt — (xx)/ m - The results are plotted using dimensionless 
quantities in Fig. f|. We see that the evidence for universal scaling behaviour 
in the dimensionless plots is much less convincing, particularly for xi f° r 
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Figure 4: /?gX« (filled) and /3(?Xt (open) vs. Pc m from a 16 3 lattice. 

$G m > 0.01. At the smallest m values the two susceptibilities even coincide 
within errors, suggesting that chiral symmetry is unbroken. This is probably 
because the data from different (3q are taken on different physical volumes, 
and two-point functions ( pTQ are prone to finite volume effects, particularly 
since a Goldstone boson is anticipated in the transverse channel as m — » 0. 
We can observe a tendency for xt to increase as this limit is approached. 
Bearing in mind the linear fit to the data of Fig. |3|, however, we might expect 
to have to go to mass values as small as Pc m ~ 0.002 before the divergence 
in Xt due to a Goldstone pole becomes dominant. Results in this regime 
in both continuum and thermodynamic limits will be required for universal 
statements about the nature of the IR limit to become feasible. 



4.2 B^0 

In this subsection we present the behaviour of fermionic matter under the 
combined action of the external magnetic field and the quantum gauge field. 
Before going on with the specific features of these results, let us remark 
that to facilitate comparison with the analytic results || we measured the 
magnetic field in units of its maximal value: thus we used the parameter 
b G [0, 1] defined by b = B/B max = B/n. 

A comparison of the dynamical and the quenched results is made in fig- 
ure |5( A 16 3 lattice has been used for the calculation, while the external 
magnetic field has been set to b = 0.094, a typical value. Measurements have 
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Figure 5: Fermion condensate versus Pa f° r quenched and dynamical 
fermions at b=0.094. 



been done at several values of the bare mass m and what we show in figure 
|5| is the extrapolation to the m — > limit. The method of extrapolation 
will be explained in more detail in the discussion of figure [| below. The 
volume dependence is not very big in the quenched case We observe 

that the values for the condensate in the quenched case are rather large for 
the relatively large value (3q = 1.0 for the gauge coupling. In principle the 
condensate should vanish for (5q —>■ oo, since the system moves to the "free" 
case, and the lattice volume is fixed. Presumably, even the quenched conden- 
sate will vanish, but the figure suggests that this will require a much bigger 
value of (3q. 

The fact that the full condensate for dynamical fermions is smaller than 
the corresponding quantity for quenched fermions has presumably to do with 
Pauli repulsion, which has clearly an effect in the dynamical fermion case. 

The next set consists of measurements of the fermion condensate versus 
the magnetic field for a 16 3 lattice at a fixed, strong coupling for the statistical 
gauge field — 0.10), for five values of the bare mass (figure^). For all five 
masses the plot consists of two parts with different behaviour. For relatively 
small b we find a dependence of the condensate on the external magnetic field 
which is nearly linear. For relatively big magnetic fields we find points that 
have a qualitatively different behaviour. The two regions merge smoothly. 
The "linear" part of the graph appears to extend in a wider interval of b 
for decreasing mass. It should be stressed that the mass dependence of the 



condensate is in sharp disagreement with the corresponding result in QT7 
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Figure 6: Fermion condensate at strong coupling versus the magnetic field 
for various bare masses. 
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where quenched fermions have been used. The condensate increases with 
the bare mass for all values of the magnetic field strength, whereas in the 
quenched case it was decreasing. The m — > extrapolation of the results lies 
probably under the lowest curve, however at this stage our statistics are not 
good enough to get reliable numbers for the limiting values of the condensate. 

It is interesting to know what happens as one moves towards weaker cou- 
plings. Based on the results of the quenched case |HJ we expect that the 
almost linear part referred to previously would be restricted to very small 
magnetic fields and eventually disappear. In figure [7| we show the results 
for the coupling (3q = 0.40. The linear part has indeed been restricted to 
the region b < 0.2. Moreover, we show the bare mass dependence along with 
the m — > limit of the condensate versus b. It is interesting to note that 
the data for each particular mass give a monotonic dependence, while the 
m — > limit starts having a local minimum at b ~ 0.5, which is character- 
istic of the free case, ie. with the statistical gauge field turned off, to which 
we now turn. The free case was studied in |J, where it was found that for 
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Figure 8: Approach to the free case in the limit of large fia and a fixed 
fermion mass m = 0.100. 

big enough b the condensate stopped being monotonic and started showing 
local maxima and minima. It developed the first local minimum at b = 0.5 
and then it had a succession of maxima and minima, up to b = 1. Moreover, 
there was a spectacular volume dependence. One expects, of course, that 
this free case will be reached for big enough (3q- We have seen this kind of 



approach to the free model in the quenched fermion case |T7| . In figure [8] we 



show the results for (5g = 0.40, 0.50 and 1.00 for two volumes, 8 and 16 
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Figure 9: Extrapolation to the zero mass limit for (3q = 0.33 and b=0.094. 
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Figure 10: /^-dependence of the condensate for two values of the magnetic 
field: b = and b = 0.094. 

The coupling (3 G = 0.50 is the point where the curve shows the first sign of 
non- monotonous behaviour. At (5q = 1.00 the successive maxima and min- 
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ima are clear. These features are shared with the quenched case. However, 
in the quenched calculation the volume dependence was not visible, while in 
this case, one can see the volume dependence at all values of the parameters. 
In the intermediate coupling regime {(3q = 0.40 and (3q = 0.50) the volume 
dependence, although clear, is rather small; at (3g = 1-00 the volume depen- 
dence is comparable with that of the free case. In view of this behaviour we 
may be confident that, at this large j3c, the limit of switching the gauge field 
off can be reached if one uses dynamical fermions, in contradistinction with 
the quenched calculation ||17|| . One should remark that in the free case the 
role of the bare mass is very important, since it is eventually the only source 
of mass generation. This is at the root of the large volume dependence show- 
ing up in the free case: at fixed volume the condensate goes over to zero for 
vanishing bare mass. In the quenched interacting model, though, the inter- 
action with the gauge field generates a dynamical mass, independently from 
the value of the bare mass. The volume dependence is thus small, permit- 
ting a smooth transition to the thermodynamic, as well as to the massless, 



limit ||23|| . In the present case with the dynamical fermions, the dynamically 
generated mass is suppressed, presumably due to the Pauli repulsion which 
is in full operation here, so the value of fto for which we have the transition 
to the free case is rather small compared to the quenched calculation. 

Now we are going to explain in some detail the procedure we followed 
to obtain the m — > limit. The simulations have been done at non-zero 
values of the bare mass; the massless limit has been taken by extrapolating 
the results for several bare masses to the limit m — > 0. Figure ||] shows the 
process of this extrapolation for fie = 0.33 and the external magnetic field b 
set equal to 0.094. The general behaviour is similar to that of the quenched 
case, in the sense that in the strong coupling region this curve is an almost 
straight line with a small slope and in the weak gauge coupling a second 
order polynomial proves necessary to perform the fit. At strong coupling 
the mass dependence is not very big, because it is the strong gauge coupling 
which dominates in the formation of the condensate. The error bars appear 
to be smaller here as compared to the gauge couplings of similar magnitude 
in the quenched case [[17} . 

Figure [H] contains the zero mass limit of the condensate (obtained through 
the procedure illustrated in figure |9|) versus (3g for b = and b = 0.094. In 
the former case the bulk of the results come from 8 3 lattice, however we also 
present some results from 16 3 , just to get some feeling about the volume de- 
pendence. We observe that in the strong coupling region the 6-dependence of 
the condensate is rather weak; on the contrary, at weak coupling, the exter- 
nal field is the main generator of the condensate, and we find an increasingly 
big 6-dependence as we move to large fie- In particular for zero external field 
the result is very small for (3g > 0.5, while the external field at these values 
of the gauge coupling still generates a relatively big condensate. 

The volume dependence of the condensate at b = 0.094 for a wide range 
of /3g is shown in figure [IT]. For the smallest volume 8 3 the data exhibit 
a rather big cusp at (3g — 0.30. For the two bigger volumes 16 3 and 20 3 
the slopes are discontinuous at /3q — 0.53 and (3g — 0.60 respectively, but 
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Figure 11: Volume dependence and finite size effects for dynamical fermions 
in the zero mass limit. 
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Figure 12: Finite size effects for dynamical fermions in the zero mass limit 
using dimensionless variables. 
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the discontinuity is smaller. It appears to be a finite size effect, possibly 
a shadow of the deconfining transition for each volume at weak coupling. 
The dynamically generated mass is smaller for the model with dynamical 
fermions as compared to the quenched model for the same value of the gauge 
coupling; so the finite size effects become more manifest. 



In figure [L^ we rescale the axes of figure [11], such that we get dimen 



sionless units. In particular, we use (3qB — > eB p h ys /g A for the x axis and 
Phixx) ~^ (xx)\phys/ 9* for the y axis. What was a cusp in figure [11] has 
been transformed into a maximum in dimensionless units. The condensate 
starts increasing with the external magnetic field, now that it is properly 
rescaled, but finally the finite size effects take over and the condensate starts 
decreasing. The value of the magnetic field where this happens increases with 
the lattice volume and presumably tends to infinity in the thermodynamic 
limit. 



5 Summary 

After reviewing and contrasting the formulation and global symmetries of 
QED 3 in continuum field theory and on the lattice, we have presented numer- 
ical results from a Monte Carlo simulation study using dynamical fermions. 
For B = our calculations used the largest system sizes studied to date to 
our knowledge; for B ^ this is the first such study. Our principal results 
are as follows: 

• For B = we have evidence, on the assumption of a linear scaling 
relation, for spontaneously broken chiral symmetry in the chiral and 
thermodynamic limits for Nf = 2 flavors of four- component fermion. 
The condensate (xx) nas an unusually small value, O(10~ 3 ) in natural 
units, which perhaps explains why this issue has proved a "difficult" 
problem in non-perturbative field theory. Our results should inform fu- 
ture studies on the system sizes needed in order for universal statements 
about the continuum theory to be possible. 

• We have studied the response of the condensate to an external magnetic 
field in the regimes of strong, intermediate and weak U(l)s coupling. 
In all cases (xx) is smaller than that found in comparable quenched 
simulations. As in the quenched case, for strong and intermediate cou- 
pling regimes the condensate increases monotonically with B, but the 
behaviour becomes non-monotonic at weaker coupling. In the weak 
coupling limit contact with the behaviour of free fermions including 
volume dependence is recovered, in contrast with the quenched model. 

• We have used our results in the chiral limit to make a plot of con- 
densate against external field in dimensionless units (Fig. [12]) , which 
enables a significant finite volume effect to be identified. The conden- 
sate increases as a smooth monotonic function of (3qB with negative 
curvature. Future studies with different values of the lattice magnetic 
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field parameter b will help to establish whether this behaviour is uni- 
versal, and hence characterises the continuum limit. We note that in 
the Schwinger-Dyson approach [p4] the regime of small magnetic fields 
could not be efficiently treated, while this is feasible with the Monte 
Carlo approach. 
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